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. Aims. To investigate the astrometric effects of stellar surface structures as a practical limitation to ultra-high-precision astrometry, 

e.g. in the context of exoplanet searches, and to quantify the expected effects in different regions of the HR-diagram. 
Methods. Stellar surface structures (spots, plages, granulation, non-radial oscillations) are likely to produce fluctuations in the in- 
tegrated flux and radial velocity of the star, as well as a variation of the observed photocentre, i.e. astrometric jitter We use the- 
■ oretical considerations supported by Monte Carlo simulations (using a starspot model) to derive statistical relations between the 

corresponding astrometric, photometric, and radial-velocity effects. Based on these relations, the more easily observed photometric 
and radial-velocity variations can be used to predict the expected size of the astrometric jitter Also the third moment of the brightness 
distribution, interferometrically observable as closure phase, contains information about the astrometric jitter 

Results. For most stellar types the astrometric jitter due to stellar surface structures is expected to be of order 10 micro- AU or greater 



o 



^ ' This is more than the astrometric displacement typically caused by an Earth-size exoplanet in the habitable zone, which is about 

' 1—4 micro-AU for long-lived main-sequence stars. Only for stars with extremely low photometric variability (< 0.5 mmag) and low 

magnetic activity, comparable to that of the Sun, will the astrometric jitter be of order 1 micro-AU, sufficient to allow the astrometric 
detection of an Earth-sized planet in the habitable zone. While stellar surface structure may thus seriously impair the astrometric 
^ ■ detection of small exoplanets, it has in general negligible impact on the detection of large (Jupiter-size) planets and on the determi- 

' nation of stellar parallax and proper motion. From the starspot model we also conclude that the commonly used spot filling factor is 

I not the most relevant parameter for quantifying the spottiness in terms of the resulting astrometric, photometric and radial-velocity 

, variations. 
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1. Introduction caused by surface structures. In some cases, e.g. for a rotating 



^, ^ . , • , spotted star, the shifts are periodic and could mimic the dy- 

The accuracy of astrometric measurements has improved ^^^^^^^^j p^jj ^ planetary companion, or even the star's par- 

tremendously in the past decades as a result of new techniques ^jj^^ ^^^-^^^ ^^^^ ^j^^^ j^^^^ ^^^^^ 

being introduced, both on the ground and in space. This devel- cuiTendy of great interest as a possible limitation of the astro- 
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H opment will continue in the next decade, e.g Gaia is to improve ^^^-^ ^^^^^^ ^^^^^^ Earth-like exoplanets. We want to 

. . ^ . parallax accuracy by another two orders of niagnitude conipared ^^^.^^^^ important these effects are for different types of 

with Hipparcos As a resuk trigonometric distances will be ob- ^^^^^^ especially in view of current an d future astro metric exo- 

tained for the Magellanic Clouds, and thousands of Jupiter-size j^^^^^ search es such as VLTI-PRIM A (iReffert et al.l l2005'). SIM 

exoplanets are likely to be found from the astrometric wobbles of pi^netQuest (IUnwinll2005h and Gaia (iLattanzi et al.l l2005). 
their parent stars. Even before that, ground-based interferometric 

techniques are expected to reach similar precisions for relative Astrometric observations determine the position of the cen- 
measurements within a small field. How far should we expect tre of gravity of the stellar light, or what we call the photocentre. 
this trend to continue? Will nanoarcsec astrometry soon be a re- This is an integrated property of the star (the first moment of the 
ality, with parallaxes measured to cosmological distances and intensity distribution across the disk), in the same sense as the to- 
Earth-size planets found wherever we look? Or will the accu- tal flux (the zeroth moment of the intensity distribution) or steUar 
racy ultimately be limited by other factors such as variable opti- spectrum (the zeroth moment as function of wavelength). In stars 
cal structure in the targets and weak microlensing in the Galactic other than the Sun, information about surface structures usually 
halo? The aim of this project is to assess the importance of such come from integrated properties such as light curves and spec- 
limitations for ultra-high-precision astrometry. In this paper we trum variations. For example, Doppler imaging (DI) has become 
consider the effects of stellar surface structures found on ordi- an established technique to map the surfaces of rapidly rotat- 
nary stars. ing, cool stars. Unfortunately, it cannot be applied to most of the 
Future high-precision astrometric observations will in many targets of interest for exoplanet searches, e.g. low-activity solar- 
cases be able to detect the very small shifts in stellar positions type stars. Optical or infrared interferometric (aperture synthe- 

sis) imaging does not have this limitation, but is with current 

Send offprint requests to: U. Eriksson, e-mail: urbanOastro .lu.se baselines (< 1 km) in practice limited to giant stars and other 
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extended objects (see Monnier et al. 2006 for a review on recent 
advances in stellar interferometry). Interferometry of marginally 
resolved stars may, however, provide some information about 
surface structures through the closure phase, which is sensitive to 
the third central moment (asymm etry) of th e stellar intensity dis - 
tribution jJVIonnie j2003l: iLachau me 2003; Labevrie et al.ll2006l) . 

Since there is limited information about surface structures 
on most types of stars, an interesting question is whether we 
can use more readily accessible photometric and spectroscopic 
data to infer something about possible astrometric effects. For 
example, dark or bright spots on a rotating star will in general 
cause periodic variations both in the integrated flux and in the 
radial velocity of the star, as well as in the photocentre and the 
asymmetry of the intensity distribution. Thus, we should at least 
expect the astrometric effect to be statistically related to the other 
effects. 

We show that there are in fact relatively well-defined statis- 
tical relations between variations in the photocentre, total flux, 
closure phase and radial velocity for a wide range of possible 
surface phenomena. These relations are in the following used to 
predict the astrometric jitter in various types of stars, without 
any detailed knowledge of their actual surface structures. 



2. Astrometric limits from previous studies 

The discovery of exoplanets by means of high-precision radial 
velocity measurements has triggered an interest in how astro- 
physical phenomena such as magnetic activity and convective 
motions might affect the observed velocities (Saaret al. 2003). 
Evidence for dark spots have been seen photometrically and 
spectroscopically for many cool stars other than the Sun, and 
quantified in terms of an empirically determined spot filling fac- 
to^, ranging from <K 1% for old, inactive stars to several 
percent for active stars. It is therefore natural to relate the ex- 
pec ted radial- velocity effect s to the spot filling factor For exam- 
ple, ISaar & Donahud (Il997h used a simple model consisting of 
a single black equatorial spot on a rotating solar-like star to de- 
rive the following relation between / (in percent), the projected 
rotational velocity V sin / and the amphtude Av, of the resulting 
radial velocity variations: 



Ay, = 0.0065 sin/ 



(1) 



In a similar vein, lHatzesl(l2002h estimated both the radial veloc- 
ity amplitude and the corresponding astrometric effect from a 
similar model, but assuming a fixed spot size (2° radius) and in- 
stead varying the number of spots placed randomly on the stellar 
surface centred around the equator. For the radial velocity am- 
plitude they found 



Ay, = (0.0086 y sin / - 0.0016) /' 



0.9 



(2) 



in approximate agreement with ([T]i, while the total amplitude of 
the astrometric effect (converted to linear distance) was 



Apos = (7.1x 10'^ AU)/ 



•0.92 



(3) 



iRefFert et al.l (|2005) discuss the accuracy and limitations of 
the PRIMA (Phase-Referenced Imaging and Micro-Arcsecond 
Astrometry) facility at the VLT Interferometer in the context of 
the search for suitable targets for exoplanetary searches, refer- 
ence and calibrations stars. According to their calculations, a 
spot filling factor of / = 2% would move the photocentre of 



a GOV star by about 3 x 10"^ AU, roughly a factor 4 less than ac- 
cording to ([3j. They also conclude that the corresponding bright- 
ness variation is less than 2%. 

But / alone may not be a very good way to quantify the 
'spottiness'. For example, the photometric or astrometric effects 
of a large single spot are obviously very different from those 
of a surface scattered with many small spots, although the spot 
filling factor may be the same in the two cases. Therefore, more 
detailed (or more general) models may be required to explore 
the plausible ranges o f the astrometric effects. 

iBastian & Hefeld (^005) give an assessment of the astro- 
metric effects of starspots, and conclude that they are hard to 
quantify, mostly because of the insufficient statistics. Although 
starspots are common among cool stars with outer convective 
zones, data are strongly biased towards very active stars. They 
conclude that the effects on solar-type stars are likely to be neg- 
ligible for Gaia, while much larger spots on K giants may be- 
come detectable. For supergiants and M giants, having radii of 
the order of lOOT?© (or more), the effect may reach 0.25 AU (or 
more), which could confuse the measurement of parallax and 
pro per motio n. 

ISozzettil (|2005) gives an interesting review of the astromet- 
ric methods to identify and characterize extrasolar planets. As 
an example of the astrophysical noise sources affecting the as- 
trometric measurements, he considers a distribution of spots on 
the surface of a pre-mainsequence (T Tauri) star. For a star with 
radius IRq seen at a distance of 140 pc, he finds that a variation 
of the flux in the visual by AF/F - 10% (rms) corresponds to 
an astrometric variation of ~ 3 yuas (rms), and that the two effects 
are roughly proportional. 

While the astrometric effects cannot yet be tested ob- 
servationally, it is possible to correlate t he photometric and 
radial-velocity variations for some stars dOueloz et al.l 120011; 
Henry et al. 2 002). From a small sample of Hyades stars 
iPaulson et al.l (l2004bl) found an approximately linear relation 



o-,,^ ^ 2 H- 3600 o-„, [ms-i] 



(4) 



' / is interpreted as the fraction of the visible hemisphere of the star 
covered by spots. 



between the RMS scatter in Stromgren y magnitude (cr,„) and 
in radial velocity (o-y^). This relation supports the idea that a 
large part of the radial- velocity scatter in these stars is caused by 
sur face structures. 

ISvensson & Ludwi3 (|2005|) have computed hydrodynami- 
cal model atmospheres for a range of stellar types, predicting 
both the photometric and astrometric jitter caused by granula- 
tion. They find that the computed astrometric jitter is almost 
entirely determined by the surface gravity g of the atmosphere 
model, and is proportional to ^ ' for a wide range of models. 
This relationship is explained by the increased granular cell size 
with increasing pressure scale height or decreasing g. The ra- 
dius of the star does not enter the relation, except via g, since 
the increased leverage of a large stellar disk is compensated by 
the averaging over more granulation cells. For their most ex- 
treme model, a bright red giant with log^ - IJRJRq ^ 95) 
they find o-p„, ^ 300 //AU. luudwig & BeckersI (l2005l) extended 
this by considering the effects of granulation on interferometric 
observations of red supergiants. They show that both visibili- 
ties and closure phases may carry clear signatures of deviations 
from circular symmetry for this type of stars, and conclude that 
convection-related surface structures may thus be observable us- 
ing interferometry. 

Ludwig ( 2006) outlines a statistical procedure to characterise 
the photometric and astrometric effects of granulation-related 
micro-variability in hydrodynamical simulations of convective 
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stars. Based on statistical assumptions similar to our model in 
Appendix A, he finds the relation 



variation of contrast) we find 



(5) 



between the RMS fluctuation of the photocentre in one coordi- 
nate (x), the radius of the star (R), and the relative fluctuations of 
the observed flux (F). 



3. Modeling astrometric displacements 

3.1. Relations for the astrometric jitter 

In a coordinate system xyz with origin at the centre of the star 
and +z away from the observer, let I{r, t) be the instantaneous 
surface brightness of the star at point r - {x, y, z) on the visi- 
ble surface, i.e. the specific intensity in the direction of the ob- 
server. We are interested in the integrated properties: total flux 
F{t), photocentre offsets Ajic(f), A3'(f) in the directions perpen- 
dicular to the line of sight, the third central moment of the in- 
tensity distribution n^it), and the radial velocity offset AvR(f). 
These are given by the following integrals over the visible sur- 
face S {z< 0): 



Fit) 
Ax(f) 
^y{t) 

^l3(t) 

AvR(f) 



r I(r,t)[x-Ax(t)ftidS 
F(t) Js 

-^T r I(r,t)[(cjxr)-z]fidS 
F(t) Js 



(6) 
(7) 
(8) 
(9) 
(10) 



where fi = \z\/R is the geometrical projection factor applied to 
the surface element when projected onto the sky, co is the angu- 
lar velocity of the star and z the unit vector along +z. (For the 
third moment, only the pure x component is considered above.) 
Equation (fTOl i assumes that the star rotates as a rigid body, that 
rotation is the only cause of the radial-velocity offset, and that 
the overall offset can be calculated as the intensity-weighted 
mean value of the local off'set across the surface. The flux varia- 
tion expressed in magnitudes is 



Am(f) = 1.086 



F(t) - (F) 



(11) 



where {F} is the time-averaged flux. 

Using a similar statistical method as iLudwigl (l2006l) . the 
RMS variations (dispersions) of m{t), Ax(t), Ay{t) and fisit) can 
be estimated from fairly general assumptions about the surface 
brightness fluctuations (Appendix |A]i. This calculation is ap- 
proximately valid whether the fluctuations are caused by dark 
or bright spots, granulation, or a combination of all three, and 
whether or not the time variation is caused by the rotation of 
the star or by the changing brightness distribution over the sur- 
face. The result is a set of proportionality relations involving the 
radius of the star R, the limb-darkening factor a, and the centre- 
to-limb variation c of the surface structure contrast [see dA.SI l 
and ( lA.lSI l for the definition of a and c]. For a - 0.6 (typical so- 
lar limb-darkening in visible light) and c = (no centre-to-limb 



0"Av = O-p, 



Q.31()Ra-,r. 



cr^, - 0.139/?V,„ 



(12) 
(13) 



where cr^ designates the dispersion of the quantity q. 

For the radial-velocity dispersion, a similar relation can be 
derived under the previously mentioned conditions of a time- 
independent, rigidly rotating star Using that (coxr) ■ z- tOxy — 
ojyX we have 



AvR(f) - WvAy(f) - w,.Ax(f) 
and 

2 22, 22 / 2, 2\2 



(14) 



(15) 



since Ax(t) and A3'(f) are statistically uncorrected according to 
Eq. iA.l\ . Noting that R(cl)1 + w^)''^ equals the projected rota- 
tional velocity V sin / we can also write ( fT5l ) as 



o-pos = Ro- „^/(V sin i) 



(16) 



which may be used to predict the astrometric jitter from the ra- 
dial velocity variations, if the latter are mainly caused by rota- 
tional modulation. Combined with ( fT2b we find under the same 
assumption 



cr,,„ - 0.376ysin/cr„ 



(17) 



In terms of the rotation period P = In/cj, and assuming random 
orientation of m in space, Eq. ( fT6] l can be written 



3 P 



= ^^^^^.-0.l95P..^ 



3.2. Modeling discrete spots 



(18) 



As a check of the general relations in Sect. 13.11 we have made 
numerical simulations with a very simple model, consisting of a 
limited number of (dark or bright) spots on the surface of a ro- 
tating star The behaviour of the integrated properties are readily 
understood in this case (cf. Fig.lTJ: 

- the flux is reduced in proportion to the total projected area of 
the visible spots (or the spot filling factor /); 

- a black spot on, say, the +x side of the star will shift the pho- 
tocentre in the -x direction and cause a negative skewness 
of the flux distribution along the x direction; 

- the apparent radial velocity of the star is modified, depending 
on whether the dark spot is located on the part of the disk 
moving towards the obse rver (g iving Avr > 0) or away from 
the observer (Avr < 0) (iGravli2005l p. 496 and references 
therein). 

Bright spots cause similar effects but with the opposite sign. 
Limb darkening of the stellar disk and a possible centre-to-limb 
variation of spot contrast will modify the precise amount of these 
shifts, but not their qualitative behaviour 

We assume a spherical star with spots that are: 

- absolutely black, 

- small compared to the stellar radius R, 

- of equal area A (measured as a fraction of the total surface), 

- randomly spread over the whole stellar surface, and 

- fixed in position on the surface, while the star rotates. 



4 



U. Eriksson and L. Lindegren: Limits of ultra-high-precision optical astrometry: Stellar surface structures 



0.015 


_l 


' 1 ' ' 


; . 1 . 1 1 




0.010 


- 






^Displacement in x 


0.005 




Radial v^ocHy 


/' 


y''^ .3rd central moment in x _ 










-''^\ 


0.000 








-0.005 


- 




\ ^-"^ ^ ■ '' 


Displacement in y _ 


-0.010 








^_,Magnitude _ 


-0.015 




1 


1 , . , 1 


1 , , , |- 




0.O 


0.2 


0.4 0.6 


0.8 1.0 



100 



Rotation Phase 

Fig. 1. The curves show the effects in magnitude, position, ra- 
dial velocity and intensity skewness (third central moment) of a 
single dark spot located at latitude 30°. The star is observed at 
inclination / = 90° and the limb-darkening parameter a = 0.6. 
The vertical scale is in arbitrary units for the different effects. 



For circular spots of angular radius p (as seen from the centre 
of the star), we have A - sin^(p/2). The assumption of abso- 
lutely black spots is uncritical if we interpret A as the equivalent 
area of the spot, i.e. the area of a completely black spot causing 
the same drop in flux. Bright spots can formally be handled by 
allowing negative A. 

The star is assumed to rotate as a rigid body with period P 
around an axis that is tilted an angle / to the line of sight (+z). 
For the present experiments we take the +y direction to coincide 
with the projection of the rotation vector cj onto the sky; thus 
tOx - Q, cjy - LJ sin and o), - lo cos where w = 27r/P. Limb 
darkening of the form intensity oc 1 - a -h OyU is assumed, where 
= \z\IR. 

To model a rotating spotted star, we place the spots of the 
given size A randomly on the surface of a spherical star and tilt 
the axis to a certain inclination /. Letting the star rotate around 
its axis we calculate the integrated quantities as functions of the 
rotational phase, taking into account the projection effect on the 
area of each spot (by the factor ju) as well the limb-darkening 
law. 

The effects of a single black spot as function of the rotational 
phase are illustrated in Fig. [T] It can be noted that the effects 
are not unrelated to each other; for example, the radial-velocity 
curve mirrors the displacement in x, and both of these curves 
look like the derivative of the photometric curve. This is not a 
coincidence but can be understood from fairly general relations 
like ( fT4l i. With many spots the curves become quite complicated, 
but some of the basic relationships between them remain. 

The total equivalent area of the spots is AN (the spot filling 
factor / ^ 2AN). As long as AN <K 1, all the effects are propor- 
tional to A. The dependence on is more complex because of 
the random distribution of spots. For example, the photometric 
effect will mainly depend on the actual number of spots k visible 
at any time. For any random realization of the model, k follows 
a binomial distribution with parameters p = 0.5 and A^; its dis- 
persion is therefore /2. We can therefore expect the RMS 
photometric effect to be roughly proportional to A y/N. Similar 
arguments (with the same result) can be made for the other ef- 
fects. 
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Fig. 2. Results of Monte Carlo simulations of rotating stars with 
different number (A^) of spots, all of the same size (A = 0.0025). 
The different graphs refer to (from top to bottom) cr„,, cTpos, cr^, 
and cr,,|j, expressed on an arbitrary scale; the dots and error bars 
show the mean value and dispersion of the cr values for a set 
of simulations with given A^. The dashed lines have slope 0.5, 
corresponding to cr oc 



Monte Carlo simulations of a large number of cases with 
A - 0.0025 (spot radius p - 5.73°) and A^ in the range from 1 to 
50 (assuming random orientation of the rotation axis and a limb- 
darkening parameter a = 0.6) indeed show that the RMS effects 
in magnitude, photocentre displacements, third central moment 
and radial velocity are all, in a statistical sense, proportional to 
VA^ (Fig.|2]i. More precisely we find 



o-,„ ^ (L17±0.60)-A Va^ 
(0.57 ± 0.25) ■ A Va^ ■ 
(0.22 ± 0.09) ■A^/N■R^ 



' pos 



cr,^ ^ (0.51 +0.26)-A Va^-/;w 



(19) 
(20) 
(21) 
(22) 



where the values after + show the RMS dispersion of the pro- 
portionality factor found among the different simulations. 

The relations (fT9]l-(l22]i suggest that a measurement of any 
one of the four dispersions can be used to statistically predict the 
other three dispersions, assuming that we know the approximate 
radius and rotation period of the star, and that the different effects 
are indeed caused by the rotating spotted surface. An important 
point is that it is not necessary to know A or A^ in order to do 
this. For example, expressing the other effects in terms of the 
photometric variation we find 



^pos 



0.497;cr,„ 

0.19 R^o-^ 
0.43 R cu o-„ 



(23) 
(24) 
(25) 



Comparing these relations with the theoretical results in (fT2b- 
( fTSl l we find that the numerical factors from the numerical ex- 
periments are systematically some 30^0% larger than accord- 
ing to the statistical theory. This discrepancy largely vanishes if 
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the models are constrained to high inchnations (/ ^ +90°). This 
suggests that the discrepancy is mainly caused by the small val- 
ues of cr„, obtained in models with small inclinations, i.e. when 
the star is seen nearly pole-on. The differences in these factors 
are in any case well within the scatter indicated in Eq. I23l425l 
which emphasizes the statistical nature of the predictions based 
e.g. on photometric variations. 

It should also be noted that there is a considerable scatter 
between the different realisations reported in Eqs. (fT9])-(l22b. 
amounting to about 50% RMS about the mean RMS effect. Thus, 
any prediction based on either (fT2b-(fT8l) or (fT9])-(l22t is only 
valid in a statistical sense, with considerable uncertainty in any 
individual case. Nevertheless, the overall agreement between the 
results of these very different models suggests that the statistical 
relations among the different effects have a fairly general valid- 
ity. The expressions for cr^,^ are the least general in this respect, 
as they obviously break down if the structures change on a time 
scale smaller than P, or if the surface structures themselves have 
velocity fields. Equations ( fTSl i and ( fT3] ) do not depend on the 
assumption that the variability is caused by the rotation. 

When modeling spotted stars, any brightenin g effect of fac- 
ulae is often disregarded (for more details see lAamm-UlvasI 
l2005h : only the darkening effect of spots is computed. For 
the Sun, the effect of faculae is known to be comparable and 
sometimes even l arger than the darkening effect of sunspots 
(Eker et alj l200l 'Chapman' 198?; 'ChapmaiT&Miyii [19861 
[chapma n et al.ll9 92; Steinegger et al. 1996). However, since the 
general relationships, e.g. in ([T2b-([T8]l, are equally valid for 
bright and dark spots (or any mixture of them), it should still be 
possible to predict the astrometric effects from the photometric 
variations. 



3.3. Comparison with previous studies and observations 

The (near-) proportionality between the observable effects and 
the spot filling factor / oc AN expressed by Eqs. ([TJ-lO is not 
supported by our spotted model, which predicts that the effects 
are proportional to A y/N. However, for small and a filUng fac- 
tor of a few percent we have rough quantitative agreement with 
these earlier results. We note that ((2]i and ^ can be combined to 
give an approximate relation similar to dTTI l. 

Equation (O derived bv lLudwigl (|2006|) is practically identi- 
cal to our ([TZt . which is not surprising as they are based on very 
similar statistical models. 

Both the theoretical result and the result from the simulation 
for the relationship between the RMS for the radial velocity and 
the RMS for the magnitude shows a distinct relatio n and this re- 
sult i s confirmed by observations in the literature dPaulson et al.l 
[2004bi) for a very limited number of stars in the Hyades all hav- 
ing rotation period of P ~ 8.5 days. These are G0V-G5V stars 
and should therefore have approximately the same radii as the 
Sun (R ~1 X 10^ km). Equation ((25]l then gives 



(T,,R ^ 2600 (T,„ [m s 



(26) 



in reasonable agreement with the empirical result in ((4|. The 
simulations by ISozzettil (|2005|) give an astrometric jitter that is 
roughly a factor 2 greater than predicted by ([T2l i or ((23). 

Thus the results of previous studies generally agree within 
a factor 2 or better with the theoretical formulae derived in this 
Section. 



4. Application to real stars 

In this section we use known statistics about the photometric and 
radial-velocity variations of real stars in order to predict the ex- 
pected astrometric jitter for different types of stars. Rather than 
using angular units, we consistently express the astrometric jitter 
in linear units, using the astronomical unit AU, mAU (lO"-' AU) 
or yuAU (10^^ AU). This eliminates the dependence on the dis- 
tance to the star, while providing simple conversion to angular 
units: 1 fiAV corresponds to 1 fias at a distance of 1 pc. We also 
note that 1 mAU - 0.215/;o and 1 fiAV ^ 150 km. 



4.1. Pre-Main Sequence (T Tauri) stars 

T Tauri stars are low-mass, pre-main sequence stars in a dy- 
namic stage of evolution often characterised by prominent dark 
spots, bipolar outflows or jets, accreting matter with associ- 
ated rapid brightness variations, and in many cases circumstellar 
disks extending to a few hundred AU (e.g., Rhode et aT]|200II : 
He rbst et al.ll2002l: ISiciha-Aguilar et al.|[2005h . Taking the star- 
forming region in the Orion nebula as an example, the spectral 
types rang e from G6 to M6, w ith the large majority in the range 
KO to M4 (iRhode et al.ll200ll ). 

Many processes may contribute to the astrometric jitter of 
these stars besides their surface structures, e.g. photometric ir- 
regularities of the circumstellar disk. The statistical relations 



derived in Sect. [3] could the refore mainly set a lower limit to 
the likely astrometric effects. iHerbst et alj (1 1 994 ) found that the 
photometric variability of (weak) T Tauri stars (WTTS) is of 
the order of 0.8 mag due to cool spots and occasion al flares. 
Assuming a typical radius of ~ 2Rq dRhode et al.l200lb . Eq. ( 1231 ) 
leads to an estimated astrometric variability of the order of 
IRe ~ 5000 yuAU. 



4.2. Main-Sequence stars 



lEver & GrenonI (1 19971) have used the Hipparcos photometric 
data to map the intrinsic variability of stars across the HR dia- 
gram. On the main sequence (luminosity class V), stars of spec- 
tral type B8-A5 and F1-F8 are among the most stable ones, 
with a mean intrinsic variability cr„, < 2 mmag and with only 
a few percent of the stars having amplitudes above 0.05 mag. 
Early B type stars are nearly all variable with a mean intrinsic 
variability of ~ 10 mmag, and among the cool stars the level 
and frequency of variability increases from late G to early M 
dwarfs. In the instability strip (A6-F0) the main-sequence stars 
are mostly micro-variable with cr„, up to several mmag. Among 
F-K stars the degree of variability is pro bably also a stron g func- 
tion of age or chromospheric activity (iFekel et al.ll2004l) : e.g., 
the Hyades (age ~ 6 00 Myr) show variations of about 10 mmag 
(iRadick et alj(l995l) . 

The Sun (G2V) is located in one of the photometrically most 
stable parts of the main sequence, and is one of the (as yet) 
few stars for which the micro-variability has been studied in de- 
tail. Analysis of the VIRGO/SoHO total solar irradiance data 
(Lanza et al. 2003) show variability at the level cr„, ^ 0.25 mmag 
(relative variance 5 x 10"^) on time scales < 30 days, which 
can largely be attributed to rotational modulation. The longer- 
term, solar-cycle related variations are of a similar magnitude. 
The optical data show a strong wavelength dependence, with 
cr„, ^ 0.2 mmag at 860 n m increasing to 0.4 mmag at 550 nm 
and 0.5 mmag at 400 nm dLanza et alj|2004 . For comparison, 
a single large sunspot group (equivalent area A = 0.05%, cor- 
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responding to / = 0.1%) gives cr,„ - 0.6 mmag according to 

03. 

The photometric variations of the Sun on short (rotation- 
related) timescales appears to be representative fo r solar-Uke 
stars of similar age and chromospheric activity jFekel et al.l 
|2004|) . Thus, we may expect cr,„ < 1 mma g for 'solar twins' can - 
didates, such as the sample studied by M elendez et alj (|2006|). 
Inspe ction of the Hipparcos photometry for these stars (lESAl 
Il997h confirm that most of them show no sign of variability at 
the sensitivity limit of a few mmag. Much more detailed and ac- 
curate statistics on micro-variability in solar-type stars are soon 
to be expected as a result of survey missions such as MOST 



([Walker et al.l l2003l) . COROT (BagUn et al. 2002) and Kepler 



Jfiasri et al.ll2005l^ 



The increased frequency and amplitude of variations for 
late G-type and cool er dwarf star s is a t least partly at- 
tributable to starspots. lAigrain et al.l (12004) estimated stellar 
micro-variability as function of age and colour index from a 
scaling of the solar irradiance power spectrum based on the 
predicted chromospheric activity level. For example, they find 
cr,„ ^1.5 mmag in white light for old (~ 4.5 Gyr) F5-K5 stars, 
practically independent of spectral type, while for young stars 
(~ 625 Myr) cr,„ increases from 2 to 7 mmag in the same spec- 
tral range. 

V ariability among field M dwarfs has been studied e.g. by 
iRoc kenfeller et al. (2006), who find that a third of the stars in 
their sample of M2-M9 dwarfs are variable at the level of cr,„ ~ 
20 mmag. Evidence for large spots has been found for many K 
and M stars, yielding brightness amplitudes of up to a few tenths 
of a magnitude. 

A large body of data on radial- velocity jitter in (mainly) F, G 
and K stars has been assembled from the several on-going planet 
search programmes and can be used to make statistical predic- 
tions as function of colour, chromospheric activity and evolu- 
tionary stage. However, since at least part of the radial-velocity 
jitter is caused by other effects than the rotation of an inho- 
mogeneous surface (e.g., by atmospheric convective motions), 
its interpretation in terms of astrometric jitter is not straight- 
forward. From the observations of ~4 50 stars in th e California 
and Carnegie Planet Search Program, Wright (2005) finds a ra- 
dial velocity jitter of ~ 4 m s ' for inactive dwarf stars of spectral 
type F5 or later, increasing t o some 10 m s~' fo r stars that are ei- 
ther active or more evolved. ISaar et al.1 (1 19981) . using data from 
the Lick planetary survey, find intrinsic radial- velocity jitters of 
2-100 m s"' depending mainly on rotational velocity (Vsin/) 
and colour, with a minimum around B-V ^ 1.0-1.3 (spec- 
tral type ~K5 ) . For a sample of Hyades F5 to M2 dwarf stars, 
iPaulson et all ( l2004al) find an average rms radial velocity jitter 
of ~16 m s"'. 



4.3. Giant stars 

For giants of luminosity class III, Hipparcos photometry has 
shown a considerable range in the typical degree of var iabil- 
ity depending on the spectral type (lEver & GrenonllT997h . The 
most stable giants (cr„, < 2 mmag) are the early A and late 
G types. The most unstable ones are of type K8 or later, with 
a steadily increasing variability up to ~ 0. 1 mag for late M gi- 
ants. The stars in the instability strip (roughly from AS to F6) 
are typically variable at the 5-20 mmag level. As these are pre- 
sumably mainly radially pulsating, the expected astrometric jit- 
ter is not necessarily higher than on either side of the insta- 
bility strip. This gen eral picture is confirmed by other studies. 
iJorissen et al.l (Il997h found that late G and early K giants are 



stable at the cr,„ < 6 mmag level; K3 and later types have an 
increasing level of micro-variability with a time scale of 5 to 
10 days, while b - y - 1.1 (^M2) marks the onset of large- 
amplitude variability (cr,„ > 10 mmag) typically on longer time 
scales (~ 100 days). From a larger and somew hat more sensitive 
survey of G and K giants, iHenrv et"aLl (l2000h found the small- 
est fraction of variables in the G6-K1 range, although even here 
some 20% show micro-variability at the 2-5 mmag level; giants 
later than K4 are all variable, half of them with cr,„ > 10 mmag. 
The onset of larg e-amplitude variabi lity coincides with the coro- 
nal dividing line (iHaisch et al.ll99l]) separating the earlier giants 
with a hot corona from the later types with cool stellar winds. 
This suggests that the variability mechanisms may be different 
on either side of the dividing line, with rotational modulation of 
active regions producing the micro-variability seen in many gi- 
ants earlier than K3 and pulsation being the main mechanism 
for the large r-amplitude variations in the later spectral types 

(iHenrv et al.ll2000.) . 

Several radial-velocity surveys of g iants (iFrinketal .1120011; 
ISetiawan et al.l2004tlHekker et al.l200 6') show increasing intrin- 
sic radial- velocity variability with B-V - 1.2, with a more or 
less abrupt change around B-V - 1 .2 (^ K3). Most bluer giants 
have cr,,^ - 20 m s"' while the redder ones often have variations 
of40-100ms-i. 



4.4. Bright giants and supergiants 

With increasing luminosity, variability becomes increasingly 
common among the bright giants and s upergiants (l uminosity 
class Il-Ia). The Hipparcos survey (E ver & GrenonI 19971) shows 
a typical intrinsic scatter of at least 10 mmag at most spectral 
types, and of course much more in the instability strip (including 
the cepheids) and among the red supergiants (including semireg- 
ular and irregular variables). Nevertheless there may be a few 'is- 
lands' in the upper part of the observational HR diagram where 
stable stars are to be found, in particular around G8II. 

It is clear that pulsation is a dominating variability mecha- 
nism for many of these objects. However, 'hotspots' and other 
deviations from circular symmetry has been observed in inter- 
ferometrical i mages of the surfaces of M supergiants and Mira 
varibles (e.g.. lTuthill et al.lll997llT999 ). possibly being the vis- 
ible manifestations of very large convection cells, pulsation- 
induced shock waves, patchy circumstellar extinction, or some 
other mechanism. Whatever the explanation for these asymme- 
tries may be, it is likely to produce both photometric and astro- 
metric variations, probably on time scales of months to years. 
Kiss et al. (2006) find evidence of a strong 1 // noise component 
in the power spectra of nearly all red supergiant semiregular and 
irregular variable stars in their sample, consistent with the pic- 
ture of irregular variability caused by large convection cells anal- 
ogous to the granulation-induced variability background seen for 
the Sun. 



4.5. Summary of expected astrometric jitter 

Table [1] summarises much of the data discussed in this Section 
for the main-sequence, giant and supergiant stars, and gives the 
corresponding estimates of the astrometric jitter (o-pos) based on 
theoretical formulae. These estimates are given in three columns 
labelled with the corresponding equation number: 

- Equation ( fT2l i is used to predict the positional jitter from the 
typical values of photometric variability in column cr,„. This 
is based on the assumption that the variability is due either 
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to (dark or bright) spots, granulation, or any other surface 
features that vary with time. Note that the temporal variation 
need not be related to stellar rotation. The resulting o-pos are 
probably realistic order-of-magnitude estimates except when 
the photometric variability is mainly caused by radial pulsa- 
tions. In such cases (e.g., for stars in the instability strip and 
red supergiants) the values given clearly represent upper lim- 
its to the real effect. 

- Equation ( fTSb is used to predict the astrometric effect from 
the radial-velocity variability in column 0-^.^. This is only 
valid if the radial velocity is rotationally modulated. Since 
pulsations, non-radial oscillations, convection and many 
other effects may cause radial-velocity variations without a 
corresponding astrometric effect, these estimates are again 
upper limits. Nevertheless, rotational modulation is im- 
portant among active (young) main-sequence stars and M 
dwarfs, and for these objects Eq. ( fTSb may provide correct 
order-of-magnitude estimates. 

- Finally we have included an estimate of the astrometric jitter 
based on the following equation 

o-po, = (300 fiAU) X 10'"'°s« (27) 

with \ogg taken from I Coxl (12000^. Equation ^ is de- 
rived from the inverse relati on to surface gravity g found by 
ISvensson^ife Ludwid (i2005h for a range of hydrodynamical 
model atmospheres. Although the authors warn that spheric- 
ity effects may render an extrapolation of this relation to su- 
pergiants very uncertain, we have applied it to all the stellar 
types in the table. Since it only includes the random effects of 
stellar granulation, it represents a lower limit to the expected 
astrometric jitter. 

If the estimates based on the photometric and radial-velocity es- 
timates are strictly considered as upper limits, the results in the 
table appear rather inconclusive. However, if the likely mech- 
anisms of the variabilities are also considered, it is possible to 
make some quantitative conclusions. For main-sequence A to M 
stars, the expected level of astrometric jitter is generally in the 
range 2-20 fj.AU probably depending mainly on the level of stel- 
lar activity; old, inactive stars should have less jitter (2-5 juAU). 
The Sun appears to be more stable than the typical old, solar-like 
star, but not by a large factor. The most stable giant stars are the 
late F to early K types, were the expected astrometric jitter is of 
order 25 fiAU. Late-type giants and supergiants have cTpos of a 
hundred to several thousand fiAU. 

5. Discussion 

5.1. Astrometric signature of exoplanets 

The possibility for an astrometric detection of a planet depends 
on the angular size of the star's wobble on the sky relative to the 
total noise of the measurements, including the astrophysically 
induced astrometric jitter discussed in the previous section. In 
linear measure, the size of the wobble is approximately given 
by the semi-major axis of the star's motion about the common 
centre of mass, or the astrometric signature 



a - 



M, 



Mp 
■ a ^ — a 



(28) 



M. -r .„p 

(cf. iLattanzi et al.ll2000 l who however express this as an angle), 
where Mp is the mass of the exoplanet, M, that of the star, and a 
the semi-major axis of the relative orbit. In all cases of interest 
here, Mp <sc M,, so that the second equality can be used. 



It is of interest to evaluate the astrometric signature for the 
already detected exoplanets. For most of them we only know 
Mp sin / from the radial-velocity curve, and we use this as a 
proxy for Mp. This somewhat underestimates the astrometric 
effect, but not by a large factor since the spectroscopic detec- 
tion method is strongly biased against systems with small sin /. 
Analysing the current (April 2007) data in the Extrasolar Planets 
Encyclopaedia (Schneider 2007) we find a median value a ^ 
1200 yuAU; the 10th and 90th percentiles are 15 and 10 000 fiAU. 

Future exoplanet searches using high-precision astrometric 
techniques may however primarily target planets with masses in 
the range from 1 to 10 Earth masses (MEanh - 3 x IO^^Mq) 
in the habitable zone of reasonably long-lived main-sequence 
stars (spectral type A5 and later, lifetime ^ 1 Gyr). For a star 
of luminosity L we may take the me an distance of t he habitable 
zone to be a ~ (L/L^y'^ AU ( Kasting et al.lll993l: iGould et al.l 
2003). In this mass range (~ 0.2-2 Mp) th e luminosity scales as 
Mf^ (based on data from lAndersen|[T99Th . so we find a oc M^-^^ 
and 



ff - (3 juAU) X 



Mp WM, 

MEarth/lMo 



1.25 



(29) 



For a planet of one Earth mass orbiting a main-sequence star, this 
quantity ranges from about 7 /iAU for an A5V star to 2.3 yuAU 
for spectral type KOV. 

iLopez et all (l2005h have argued that life will have time to de- 
velop also in the environments of subgiant and giant stars, during 
their slow phases of development. The habitable zone may ex- 
tend out to 22 AU for a 1 Mq star, with a correspondingly larger 
astrometric signature. However, the long period of such planets 
would make their detection difficult for other reasons. 



5.2. Exoplanet detection 

The detection probability is in reality a complicated function of 
many factors such as the number of observations, their tempo- 
ral distribution, the period and eccentricity of the orbit, and the 
adopted detection threshold (or probability of false detection). A 
very simplistic assumption might be that detection is only pos- 
sible if the RMS perturbation from the planet exceeds the RMS 
noise from other causes. Neglecting orbital eccentricity and as- 
suming that the orbital plane is randomly oriented in space, so 
that (sin^ /) = 2/3, the RMS positional excursion of the star in 
a given direction on the sky is a/ V3. With a sufficiently pow- 
erful instrument, so that other error sources can be neglected, 
the condition for detection then becomes a/o-p^^ ^ V3. In real- 
ity, a somewhat larger ratio than V3 is probably required for a 
reliable detection, especially if the period is unknown. For ex- 
ample, Sozzetti (2005) reports numerical simulations showing 
that a/cr > 2 is required for detection of planetary signatures by 
SIM or Gaia, where cr is the single-epoch measurement error, 
provided that the orbital period is less than the mission length. 
(For the corresponding problem of detecting a periodic signal in 
radial-velocity data, Marcv et al. (2005) note that a velocity pre- 
cision of 3 m s"' limits the detected velocity semi-amplitudes to 
greater than ~10 m s"', implying an even higher amplitude/noise 
ratio of 3.3.) As a rule-of-thumb, we assume that detection by the 
astrometric method is at least in principle possible if 



cr„os < 0.5q' 



(30) 



For old, solar-type stars the expected astrometric jitter is 
:^5 /iAU, implying that exoplanets around these stars with a ^ 
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Table 1. A summary of typical photometric and spectroscopic variability for different stellar types, and inferred levels of astrometric 
jitter (cTpos). The jitter is estimated in three different ways: from the photometric variability, using Eq. (fTSl i [this will overestimate 
the jitter if part of the variability is due to radial pulsation]; from the radial-velocity variability, using Eq. ( fTSb [this method will 
overestimate the jitter if the variability is not caused by rotational modulation]; and from the surface gravity, using Eq. (|27] l [this 
only includes jitter caused by granulation, and is the refore a low er limit]. References to typical observed quantities are given as 
footnotes. Radii and \ogg (not shown) are taken from lCoxl (|2000|) . 
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References: lAigrain et ai]f2QQ4h. ICoj] f2Q03). jEver & Grenor] {\99%, 'fpekel et alj flOOl. iFrink et all f20QTl. ^Hekker et 



iLanza et aTHSOOl. ^Paulson et al1l2004tD. iRadick et alj|T993). IRockenfelleretai] 12003). ISaar et al. 



120061) , iHeiirv et alj t200l. ' Ijorissen et alj|1997ft . 
T99I. 'ISetiawanetai]t2004l) 



10 //AU could generally be detected and measured astrometri- 
cally. This applies to the vast majority (>90%) of the exoplanets 
already detected by the radial-velocity method. Such observa- 
tions would be highly interesting for obtaining independent in- 
formation about these systems, in particular orbital inclinations 
and unambiguous determination of planetary masses. 

Exoplanets of about 10 MEarth orbiting old F-K main- 
sequence stars in the habitable zone {a ^ 20-50 /iAU) would 
generally be astrometrically detectable. This would also be the 
case for Earth-sized planets in similar environments {a ^ 2- 
5 yuAU), but only around stars that are unusually stable, such as 
the Sun. 

5.3. Determination of parailax and proper motion 

The primary objective of high-precision astrometric measure- 
ments, apart from exoplanet detection, is the determination of 
stellar parallax and proper motion. We consider here only briefly 
the possible effects of stellar surface structures on the determi- 
nation of these quantities. 

Stellar parallax causes an apparent motion of the star, known 
as the parallax ellipse, which is an inverted image the Earth's 
orbit as viewed from the star. The linear amplitude of the par- 
allax effect is therefore very close to 1 AU. (For a space obser- 
vatory at the Sun-Earth Lagrangian point L2, such as Gaia, the 
mean amplitude is 1.01 AU.) Thus, the size of the astrometric 
jitter expressed in AU can directly be used to estimate the min- 



imum achievable relative error in parallax. For main-sequence 
stars this relative error is less than lO"'*, for giant stars it is of 
order lO"** to lO"-', and for supergiants it may in some cases ex- 
ceed 1%. We note that a 1% relative error in parallax gives a 2% 
(0.02 mag) error in luminosity or absolute magnitude. 

If proper motions are calculated from positional data sep- 
arated by T years, the random error caused by the astrometric 
jitter, converted to transverse velocity, is ^ cTpos V2/r. Even for 
a very short temporal baseline such as T = 1 yr, this error is usu- 
ally very small: -0.1 m s"' for main-sequence stars and -0.5- 
5 ms"' for giants. (Note that 1 AU yr"' ^ 4.74 km s"' .) In most 
applications of stellar proper motions this is completely negligi- 
ble. 



6. Conclusions 

For most instruments on ground or in space, stars are still unre- 
solved or marginally resolved objects that can only be observed 
by their disk-integrated properties. The total flux, astrometric po- 
sition, effective radial velocity and closure phase are examples of 
such integrated properties. Stellar surface structures influence all 
of them in different ways. Our main conclusions are: 

1. Theoretical considerations allow to establish statistical re- 
lations between the different integrated properties of stars. 
Under certain assumptions these relations can be used to pre- 
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diet the astrometrie jitter from observed variations in pho- 
tometry, radial velocity or closure phase. 

2. The total flux, astrometrie position and third central mo- 
ments (related to closure phase) are simple moments of the 
intensity distribution over the disk, and for these the statis- 
tical relations are valid under fairly general conditions - for 
example, they hold irrespective of whether the variations are 
caused by spots on a rotating star or by the temporal evo- 
lution of granulation. By contrast, radial-velocity variations 
can only be coupled to photometric and astrometrie varia- 
tions if they are primarily caused by rotational modulation. 

3. The theoretical relations are supported by numerical simula- 
tions using a model of a rotating spotted star. In this case the 
variations in total flux, position, radial velocity and closure 
phase are all proportional to A y/N, where A is the equiva- 
lent area of each spot and the number of spots. This means 
that, e.g., the astrometrie jitter can be (statistically) predicted 
from the photometric variability without knowing A and A^. 
It is noted that the spot filling factor, being proportional to 
AN, is not the most relevant characteristic of spottiness for 
these effects. 

4. Using typical values for the observed photometric and radial- 
velocity variations in ordinary stars, we have estimated the 
expected size of the astrometrie jitter caused by surface 
structures (Table[T]l. The estimates range from below 1 fj.AU 
for the Sun, several //AU for most main-sequence stars, some 
tens of yuAU for giants, and up to several mAU for some su- 
pergiants. 

5. The expected positional jitter has implications for the possi- 
ble astrometrie detection of exoplanets. While planets heav- 
ier than 10 Earth masses may be astrometrically detected in 
the habitable zone around ordinary main-sequence stars, it is 
likely that Earth-sized planets can only be detected around 
stars that are unusually stable for their type, similar to our 
Sun. 

6. Stellar surface structures in general have negligible impact 
on other astrometrie applications, such as the determination 
of parallax and proper motion. A possible exception are su- 
pergiants, where very large and slowly varying spots or con- 
vection cells could limit the relative accuracy of parallax de- 
terminations to a few per cent. 
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Appendix A: Statistical properties of the spatial 
moments of the intensity distribution across a 
stellar disk 

In this Appendix we derive the mean values and variances of the 
moments {;c"y) for a spherical star, where x and y are spatial 
coordinates normal to the line-of-sight and () denotes the instan- 
taneous flu x-weig ht ed me an. The analysis extends and general- 
izes that of iLudwiS ( |2006|) by considering also the third moment 
(relevant for measurement of closure phase) and a centre-to-limb 
variation of the intensity contrast. 

Let 6, ip be polar coordinates on the stellar surface with -Q 
at the centre of the visible disk and ip -Q along the x axis. With 
H = cos we write the instantaneous intensity across the visible 
stellar surface S as I{fj.,(p), and introduce the non-normalized 
spatial moments 

M„„ = ^d5/(/i,^)//x'"/ = 

^m+„+2r r d^/(^,,^)^(l -ju2)('n+")/2cos"V sin^ (A.l) 

Jo Jo 

where x = R{1 - fi^y^^ costp and y = R(l - ju^)'^^ sintp. The 
factor fi in the integrand is the foreshortening of the surface ele- 
ment dS = R^dfid(p when projected normal to the line-of-sight. 
The normalized moments are given by 



M„, 



Moo 



(A.2) 



where it can be noted that Moo equals the instantaneous total 
stellar flux. 

It is assumed that I(n,tp) varies randomly both across the 
stellar surface (at a given instant), and as a function of time. As a 
consequence, the spatial moments ( lA.ll i and iA.2\ are also ran- 
dom functions of time, and the goal is to characterize them in 
terms of their mean values and variancesQ Since we are inter- 
ested in quite small effects of the surface structure it is generally 
true that the dispersions are small compared with the total flux 
and scale of the star, so that for example D[Moo] <K E[Moo] 
and D[Mi()] <K R E[M()o]. In this case the variabiHty of {x"'y"} is 
mainly produced by the numerator in ( |A.2t , and we may use the 
approximations 



E[{^y"}] = 



E[M„„] 
E[Moo] 



D[<x-/)] 



E[M, 



(A.3) 



00 J 



In the following we therefore focus on deriving the mean values 
and dispersions of the non-normalized moments M„,„ . The (tem- 
poral) mean value and dispersion of I{fi, (f) are assumed to be 
independent of ip; thus 

E[I(,x,<p)]^A(pt), m(f^,<f)]^(r,(fi) (A.4) 
where A(ju) and cr/(yu) are functions to be specified. 

A.1. Mean value of the moments 

We assume a linear limb-darkening law with coefficient a, such 
that 



A(ji) - (I - a + an)A 



(A.5) 



^ We use the notation E[X] for the mean value (expectation) of the 
generic random variable X, Y[X] = E[{X - ElX] f] for the variance and 



where Ai is the mean intensity at the disk centre (// = 1). From 
dA.ll) we obtain 



E[Af;^;j] 

R 



I r-ln 

m+n+2 I A., f A,^ A t . .\ ,.2\(m+«)/2 m 



1H 

Jo Jo 



d>pA(^i)ti{i-ny 



which with ( IA.5I) evaluates to 
E[M„„] = 

if either m or n is odd, and to 

(m- !)!!(«- 1)!! 



cos'"^' sin"^ (A.6) 



(A.7) 



E[M,„„] = 2nR 



m+n+2 



A 1 //,„+„(«) 



(m -H n) ! ! 

if both m and n are evenQ Here, we introduced the functions 

Hk(a) = f dfi(l -a + afi)fi(l -i/f'^ 
Jo 



(A.8) 



■ a (even k) 



1 — a 



k + 2 ik + 3)1 

l-a k\\ n 

+ — — TTTT — a (odd k) 



(A.9) 



(A. 10) 



(A.ll) 



(A. 12) 



I k + 2 (k + 3)\\ 2 

of which, presently, we only need 

1 fl \ la 

Hoia) and Hoia) 

^ 2 6 2^ 4 60 

Thus, the mean total flux is 
E[Moo\^nR^A, (l - ^) 
and the second moments 
E[M2o] = E[Mo2] = \nR^M (l - 
The rms extension of the star in either coordinate is given by 

\e[Moo]/ \E[Moo]/ 2 \ l-fl/3 / 
A.2. Dispersion of the moments 

In order to compute the dispersion of M„,„ we n eed to introduce 
the second-order statistics of /(yu, p). Following iLudwij (l2006h 
we divide the visible hemisphere into equal surface patches of 
size AA = /J^AeAi/) sin - R^A/jAip - 2nR^IN, with the centi-e 
of patch k at position (yu^, (pk). Thus the integral over the visible 
surface of any function gip., tp) can in the limit of large be 
replaced by a sum: 

r' r^" 2nR^ ^ 

R^\ d^l\ d^giii, ^) - — - y g{^Jiu, n) (a.14) 

Jo Jo k~l 



In particular, from (lA.lb we have 

ym+n+2 ^ 



2nR" 



N 



^ hfik (1 - ;u2)('"+«)/2 cos"-^, sinVi (A.15) 



k=i 



D[X] = V[X]''- for the rms dispersion. 



^ The double factorial notation means k\\ = k(k - 2){k - 4) • • • 2 for 
even integer k, and k\ \ = k(k - 2)ik - 4) • • ■ 1 for odd k. We have 0! ! = 
(-1)!! = 1. 
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where 4 is the mean value of if) in patch k. This expresses 
the moment as a linear combination of the random variables 4. 
If we now assume that the intensity variations A4 - h- E[//t] of 
the patches are uncorrelated, i.e. E[A/iA4'] = for k + k' , we 
have 

V[M,„„] = 



4jj.2^2»!+2«+4 JL 



N2 



y V[4] f4 (1 - fij)'"^" cos^'V, sin^V, (A.16) 



For sufficiently large the patches would resolve even the 
smallest surface structures and we would have V[4] - cr^iiUk) 
according to ( IA.4b . However, in that case the intensities of ad- 
jacent patches would be correlated, so (IA.16b would not hold. 
For the latter equation we effectively need patches that are larger 
than the correlation length of the surface structures. We must 
therefore assume that A^ is large enough for the discretization 
( |A.12t to hold, and still small enough that the patches are un- 
correlated. In this regime we have V[4] < cr^(jUjt), since 4 is 
the average intensity in patch k, not the local intensity at point 
(fik, fk)- In fact, V[4] will depend on the patch size (or A^) in 
such a way that V[4]/A^ is invariant ('Ludwig''2006). (This is 
obviously the case for independent patches: grouping them into 
larger and fewer patches decreases the variance in proportion to 
the resulting A^.) We write the invariant quantity as 

V[4] 9 9 

-Lli ^ A(fik)^C(Mkf (A. 17) 

where A(j/) is the mean intensity as before and C(ju) the centre- 
to-limb variation of the contrast (scaled by N^^^^). In analogy 
with ( IA.5b we assume a linear centre-to-limb variation of the 
contrast according to 

C(ji)^(l-c + ciu)Ci (A.18) 

Inserting (IA.17b and (lA.lSb into ( IA.16b and using (IA.14I) to 
transform the sum into an integral gives 
V[M„,„] = 27r«2'«+2«+4^ 

pi pin 

X \ djj \ d<pAQifC(fif i/ (I- //^^"cos^'" if sin^" ifi 
Jo Jo 

^ 4,2^2,,.2„.4(2^^i_^^ (^.19) 

{2m + 2n)\\ 
where we introduced the functions 

Kk{a,c)^ { An{l-a + a^lf{\-c + cnf^?(\-^l^f (A.20) 
Jo 

For A; = 0, ... 3 we have 

11 1 , , 
Koia, c) (fl -H c) H (a + 4ac + c ) 



3 6 



30 



1 



1 



2„2 



ac(a + c) H a c 

30 105 

2 1 1 , , 

11 . ^ 1 2 2 
ac(a + c) H a c 

420 126 



8 29 



23 



K2{a, c) (fl -H c) H {a + 4ac + c ) 



K^(a, c) 



105 420 
3_ 

~140 

16 13_ 

315 ~ 252 
25 



1260 
47 



2„2 



ac(a + c) -\ a c 

^ ' 6930 

29 

(a + c) + (a^ + 4ac + c^) 



1980 



38 



2 2 



^3g^.ac(a + c) + ^«c- 



(A.21) 



(A.22) 



(A.23) 



(A.24) 



Using ( IA.3b we obtain the following general expression for the 
dispersion of the normalized spatial moment D,„„ = D[{x'"y"}]: 



(2ot- 1)!! (2n- 1)! 



(2m + 2n)\ 



K,„+„(a, c) 



Ho(a) 



(A.25) 



Note that Doo - D[Moo]/E[Moo] is the relative dispersion of the 
total flux, Dio is the dispersion of the photocentre along the x 
axis, etc. We have in particular 



Dqo - C\ 



V Ki,){a, c) 



Ho(a) 



Dqi ^CiR 



D20 - D02 - C\R 



D 



30 



£>o 



C, R' 



^jjKi {a,c) 

Ho(a) 
^ lK2ia, c) 

Ho(a) 
^■^K3(a,c) 

ma) 



(A.26) 



(A.27) 



(A.28) 



(A.29) 



A.3. The third central moment 

Closure phase is sensitive to the asymmetry of the stellar image, 
and the third moments (M„„, for m + n = 3) are intended to 
provide a statistical characterization of this asymmetry. These 
moments are calculated with respect to the geometrical centre of 
the disk {a.tx-y- 0). However, intrinsic image properties such 
as size, shape and asymmetry are more properly expressed with 
respect to the photocentre, at xq - M\qIMqq, jq - Mq\ /Mqo, i.e., 
by means of central moments (here denoted with a prime). For 
example, the third central moment along the x axis is given by 



M' 



30 



I 



dS I{ii,ip)ii{x- X()f 



M30 - 3x0^20 + 3xlMio ■ 



= M30 - 3 



+ 2 



00 



■"^00 



(A.30) 



It is seen that E[M3g] = as expected. However, to calculate 
the variance of M'-^q it is necessary to make some approxima- 
tions. First we replace the even moments in the right-hand side 
of ( IA.26b by their mean values and introduce the rms extent s of 
the stellar disk from ( |A.13b . yielding 



ML = M30 - 3Mio 



(A.31) 



The photocentre displacement is normally very small compared 
with the size of the disk, so that the second term in the parenthe- 
ses can be neglected. Then 

= V[M3o] 



■ 6s^V[M2o] + 9/V[Mio] 



(A.32) 



using that E[MioM3o] - V[M2o]. Finally, the dispersion of the 
normalized third central moment is found to be 



D' 



30 



CiR^ 
~tUa) 



5 9 H2(a) 

— A3(fl, c) K^ia, c) ■ 

16 8 Ha(a) 



9 lH2(a) 



8 \Ho{a) 



K\ {a, c) 



1/2 



(A.33) 
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A.4. Scaling relations 



Since all the dispersions in ( IA.26I ) and ( IA.33I) are proportional to 
Ci we obtain a set of simple scaling relations that can be used to 
predict the dispersion of a certain moment from a measurement 
of the dispersion of a different moment (assuming that a and c 
are approximately known). The most useful relations allow to 
predict the dispersions of the first and third moments from that 
of the total flux (zeroth moment). For the first moment (photo- 
centre) we have 



D[Mio] 
D[Moo] 



R 



Ki(a,c) 



2Ko{a,c) 



1/2 



(A.34) 



and for the third central moment 



D[M^o] 
D[Moo] 



5 K-i{a,c) 9 H2{d) K2(a,c) 



16 KQ{a,c) 
9 lH2{a) 



8 Haia) Ko(a, c) 

^ 1/2 

Ki (a, c) 



8 \ //()(«)/ KQ(a,c) 



(A.35) 



The numerical factors on the right-hand sides of ( IA.34I ) and 
( IA.35I I are graphically shown in Figs. IA.ltiA.2l as functions of 
the structural parameters a and c. 




Fig.A.2. Similar to Fig. lA.ll but for the scaling factor 
D[M']/R^D[Moo] fromEq. dAJSl) 




Fig.A.l. The scahng factor D[Mn)] /RD[M oo] from Eq. ( IA34| | 
between the expected dispersions in photocentre position and to- 
tal stellar flux, plotted as function of the limb-darkening param- 
eter a and the parameter for the centre-to-limb variation of sur- 
face structure contrast (c). The different curves represent, from 
top to bottom, c=-l, -0.5, 0, 0.5, 1. The solar symbol indicates 
the typical value for solar granulation in white light, a = 0.6 and 
c = 0.4. 



